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This  report  describes  the  results  of  a  continuing  investigation  of  the  characteristics 
of  the  low  frequency  ground  motions  produced  by  atmospheric  explosions.  The  studies  re¬ 
ported  here  have  focused  on  two  aspects  of  this  problem:  (1)  the  full-scale  validation  of 
the  adopted  model  through  comparisons  with  ground  motion  data  recorded  from  atmospheric 
nuclear  tests,  and  (2)  the  preliminary  evaluation  of  a  theoretically-based  prediction 
methodology. 


With  regard  to  the  validation  study,  the  mathematical  model  has  been  applied  to  the 
theoretical  simulation  of  ground  motion  data  recorded  from  the  Nevada  Test  Site  atmospheric 
nuclear  explosions  TUMBLER  I,  UPSHOT-KNOTHOLE  10  and  SMALL  BOY.  The  results  of  these  simu¬ 
lations  have  confirmed  the  fact  that  the  observed  low  frequency  ground  motions  from  such 
explosions  can  be  accounted  for  by  the  airblast-induced ,  elastic  surface  waves.  Moreover, 


20  distribution  /  availability  of  abstract 
□  uNCLASSlf'EOAJNUMITEO  SI  SAME  AS  RPT 


□  OTIC  USERS 


21 


ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22*  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Bettv  L.  Fox _ 


22b  TELEPHONE  (ICKhidt  Are*  Code) 

(202)325-7042 _ 


22c.  OFFICE  SYMBOL 
ONA/STTI _ 


00  FORM  1473.  a*  MAR 


83  APR  edition  may  bd  used  until  etheustdd. 
All  ottidf  editions  art  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


UNCLASSIFIED 


iicumry  CLAMiricATioN  or  r»i»  »*a« 


19.  ABSTRACT  (Continued) 

the  simulation  analysis  of  the  ground  motion  data  recorded  from  the  near-surface  SMALL  BOY 
explosion  has  confirmed  the  existence  of  an  "exclusion  radius"  for  such  events,  within  which 
incident  airblast  energy  is  not  efficiently  coupled  into  the  surface  wave  mode  of  propaga¬ 
tion  due  to  dissipation  associated  with  strong,  nonlinear  interaction  effects.  On  the 
other  hand,  the  evidence  provided  by  the  TUMBLER  I  and  UPSHOT- KNOTHOLE  10  analyses  has  in¬ 
dicated  that  this  exclusion  radius  is  effectively  zero  for  explosions  with  significant 
height  of  burst. 

In  the  prediction  element  of  the  study,  the  experience  gained  with  the  Pre-DIRECT 
COURSE  experiment  has  been  assessed  and  some  approximate  scaling  laws  describing  the  de¬ 
pendence  of  low  frequency  ground  motion  characteristics  on  variables  such  as  yield,  distance, 
site  geology  and  instrument  depth  have  been  derived  and  discussed.  It  was  shown  that  pre¬ 
dicted  surface  wave  characteristics  are  quite  sensitive  to  the  near-surface  shear  wave 
velocity  structure  at  the  site  and  that  Pre-DIRECT  COURSE  experience  indicates  that  this 
velocity  distribution  cannot  be  adequately  constrained  on  the  basis  of  P  wave  refraction 
surveys.  For  a  given  site,  and  for  explosions  which  are  large  enough  that  the  surface  wave 
motions  are  dominated  by  frequency  components  near  the  corner  frequency  of  the  Rayleigh 
wave  site  response  function,  it  has  been  shown  that  the  theoretical  model  predicts  that 
the  amplitude  of  the  low  frequency  ground  motion  should  vary  with  yield,  w,  as  W1/3  and 
with  range,  R,  approximately  as  R-5/6. 
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SECTION  1 


INTRODUCTION 

As  a  result  of  research  conducted  in  recent  years,  it  is  now 
well  documented  that  the  ground  motion  displacements  observed  from  surface 
and  atmospheric  explosions  outside  the  strong  shock  regime  (i.e.,  at  ranges 
where  the  peak  incident  overpressure  has  dropped  below  about  1000  kPa)  are 
dominated  by  the  low  frequency,  surface  wave  components  of  the  motion 
(Murphy,  1978;  Murphy  and  Bennett,  1980;  Murphy  et  al_. ,  1981;  Murphy  et  al . , 
1982).  It  follows  that  in  applications  in  which  strategic  system  vulner¬ 
ability  at  such  ranges  is  an  issue  of  concern,  it  is  important  that  a 
reliable  prediction  capability  be  available  for  estimating  the  character¬ 
istics  of  these  low  frequency  ground  motions.  Significant  progress  toward 
the  development  of  such  a  capability  has  been  achieved  over  the  past  sev¬ 
eral  years,  during  which  time  a  theoretical  model  has  been  formulated  and 
successfully  applied  to  the  deterministic  simulation  of  the  low  frequency 
ground  motions  observed  from  a  variety  of  HE  tests  (e.g.,  Murphy,  1978; 
Murphy  et  al_. ,  1981).  However,  some  questions  still  remain  concerning 
the  extrapolation  of  this  HE  experience  to  the  full  scale  nuclear  environ¬ 
ments  which  are  of  interest  in  ground  motion  assessment  studies. 

The  objectives  of  the  analyses  described  in  this  report  have 
been  to  conduct  a  full  scale  validation  study  of  the  selected  model  using 
data  recorded  from  atmospheric  nuclear  explosions  at  the  Nevada  Test  Site 
(NTS)  and  to  evaluate, in  a  preliminary  fashion,  the  feasibility  of  employ¬ 
ing  a  theoretically-based  prediction  methodology  for  routine  applications. 
The  results  of  the  full  scale  validation  study  are  presented  in  Section  II, 
where  the  ground  motions  observed  from  three  atmospheric  nuclear  explo¬ 
sions  representing  a  wide  range  in  scaled  height  of  burst  (HOB)  are  theo¬ 
retically  simulated  by  computing  the  surface  wave  response  of  models  of 
the  subsurface  site  geology  to  airblast  loads  predicted  by  the  Speicher 
and  Brode  (1981)  analytic  approximations. 

Some  preliminary  considerations  concerning  the  prediction  of 
low  frequency  ground  motion  parameters  are  addressed  in  Section  III,  where 
the  experience  gained  with  the  Pre-DIRECT  COURSE  experiment  is  assessed 
and  some  approximate  scaling  laws  describing  the  dependence  of  low 


frequency  ground  motion  characteristics  on  variables  such  as  yield,  dis¬ 
tance,  site  geology  and  instrument  depth  are  presented  and  discussed. 

This  is  followed  in  Section  IV  by  a  summary,  together  with  a  statement  of 
conclusions  regarding  our  current  state  of  knowledge  concerning  the  low 
frequency  ground  motions  produced  by  atmospheric  explosions. 


SECTION  2 


ANALYSIS  OF  ATMOSPHERIC  NUCLEAR  EXPLOSION  DATA: 

FULL  SCALE  VALIDATION 

2.1  MODEL  DEFINITION. 

The  model  which  has  been  developed  to  simulate  the  low  fre¬ 
quency  ground  motions  observed  from  atmospheric  explosions  is  based  on 
the  analytic  solution  for  the  elastic  response  of  a  multilayered  model  of 
the  subsurface  site  geology  to  a  propagating  normal  load  acting  on  ^,ie 
surface  (i.e.,  an  approximation  to  an  airblast).  Assuming  that  the  air- 
blast  load  is  axisymmetric,  the  Fourier  transform  of  the  Rayleigh  wave 
component  of  the  displacement  for  a  particular  mode  at  observation  point 
r,  W(oj,r),  can  be  written  as  a  spatial  integration  over  the  Fourier  trans¬ 
form  of  the  overpressure,  p(r,a>)  in  the  form  (Murphy  and  Bennett,  1980): 

W(co,r)  =  -iirAR(<»)  H^(kr)  J*  p(r',cu)  JQ(kr')  r'  dr' 

ro  (1) 

00 

+  J0(kr)  j*  p(r",u)  H^(kr')  r'  dr' 
r 

( 2 ) 

where  AR(w)  is  the  Rayleigh  wave  site  response  function,  JQ,  Hg  denote 
the  Bessel  and  Hankel  functions  respectively,  k  is  the  wavenumber  of  the 
Rayleigh  wave  component  with  frequency  w  and  rg  denotes  the  exclusion 
radius  which  defines  the  region  within  which  strong  nonlinear  effects 
occur.  For  surface  bursts,  it  has  been  empirically  determined  that  rg 
coincides  approximately  with  the  crater  radius.  However,  for  atmospheric 
explosions  with  significant  HOB  (i.e.,  greater  than  rg),  this  heuristic 
approximation  can  be  avoided  and  rg  is  set  to  zero.  For  purposes  of  the 
present  investigation,  the  airblast  loading  function,  p(r,t),  has  been 
specified  using  Brode's  recently  revised  approximation  (Speicher  and  Brode, 
1981)  which  expresses  the  overpressure  for  a  given  yield  (W)  and  HOB  as  a 
function  of  range  and  time  through  products  of  functions  involving  fairly 
elaborate  polynomial  fits  to  the  characteristics  of  airblasts  observed 


from  surface  and  atmospheric  nuclear  explosions.  The  numerical  procedures 
required  to  transform  these  time  domain  descriptions  into  the  correspond¬ 
ing  frequency  domain  representations  required  in  the  evaluation  of  Equa¬ 
tion  (1)  (i.e.,  p(r,tj))  have  been  described  in  detail  by  Murphy  et  aK 
(1982).  Thus,  given  a  multilayered  approximation  to  the  subsurface  site 
geology  and  associated  elastic  parameters  for  each  layer,  the  airblast- 
induced  Rayleigh  wave  motions  expected  from  atmospheric  explosions  of 
arbitrary  yield  and  HOB  can  be  calculated  using  equation  (1). 

2.2  SIMULATION  ANALYSIS  OF  SELECTED  ATMOSPHERIC  NUCLEAR  EXPLOSIONS 

AT  NTS. 

The  analytical  model  described  above  has  been  successfully 
tested  against  data  recorded  from  a  variety  of  HE  simulation  tests  (e.g., 
Murphy,  1978;  Murphy  et  aK ,  1981).  However,  some  questions  still  remain 
concerning  the  extrapolation  of  this  HE  experience  to  the  full  scale 
nuclear  environments  which  are  of  interest  in  ground  motion  assessment 
studies.  For  this  reason,  a  preliminary  simulation  analysis  has  been 
conducted  in  an  attempt  to  test  the  analytical  model  against  ground  motion 
data  measured  from  NTS  atmospheric  nuclear  explosions.  The  events  which 
have  been  selected  for  analysis  are  the  TUMBLER  I,  UPSHOT-KNOTHOLE  10  and 
SMALL  BOY  explosions,  all  of  which  were  detonated  over  the  Frenchman  Flat 
testing  area  at  NTS.  The  characteristics  of  these  explosions  are  summa¬ 
rized  in  Table  1  where  it  can  be  seen  that  they  cover  a  wide  range  of 
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scaled  HOB  extending  from  242  m/kt  to  essentially  zero.  Thus,  these 
tests  provide  an  opportunity  to  assess  simulation  capability  over  essen¬ 
tially  the  entire  range  of  potential  interest  in  ground  motion  vulnera- 
bi 1 ity  studies . 

The  subsurface  geological  model  of  Frenchman  Flat  which  has 
been  adopted  for  initial  simulation  purposes  is  shown  in  Figure  1  (Sauer, 
1958;  Lipner,  personal  communication,  1981),  where  it  can  be  seen  that  it 
is  relatively  simple,  consisting  of  two  layers  over  a  halfspace.  In  this 
figure,  a  and  6  denote  the  compressional  and  shear  wave  velocities  and  p 
the  density.  The  fundamental  mode  Rayleigh  wave  dispersion  curves  corre¬ 
sponding  to  this  site  model  are  shown  in  Figure  2,  and  it  can  be  seen  that 
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there  exists  a  very  pronounced  group  velocity  minimum  at  a  frequency  of 
about  1  Hz.  Consequently,  as  is  illustrated  in  Figure  3,  the  Rayleigh 
wave  site  response  function  (i.e.,  AR(u>)  in  Equation  (1))  drops  off  sharply 
for  frequencies  below  1  Hz.  It  follows  that  for  explosions  in  the  yield 
range  of  interest  in  the  present  investigation,  it  can  be  expected  that 
the  simulated  low  frequency  ground  motions  wii1  be  dominated  by  components 
with  frequencies  of  about  1  Hz  (Murphy  and  Bennett,  1980). 

2.2.1  TUMBLER  I  . 

Selected  low  frequency  ground  motions  observed  from  this  test 
have  been  simulated  by  analytically  loading  the  surface  of  the  geologic 
model  of  Figure  1  with  an  airblast  loading  function  computed  from  the 
Speicher  and  Brode  (1981)  analytic  approximation  for  a  yield  of  1  kt  and 
an  HOB  of  242  m.  In  this  case,  the  HOB  is  sufficiently  large  that  no 
significant  nonlinear  effects  are  anticipated  and  rQ  in  Equation  (1)  is 
assumed  to  be  zero.  The  TUMBLER  I  vertical  particle  velocity  waveforms 
observed  at  ranges  of  407  m,  482  m  and  558  m  are  shown  in  Figure  4  where 
they  are  compared  with  the  corresponding  theoretical  Rayleigh  wave  syn¬ 
thetic  waveforms.  It  can  be  seen  that  the  synthetics  agree  remarkably 
well  with  the  observed  waveforms  up  to  the  point  where  the  observed  data 
have  been  arbitrarily  terminated  during  the  data  reduction.  In  particu¬ 
lar,  the  dominant  frequency  content  and  amplitude  level  of  the  major  low 
frequency  pulse  is  reproduced  to  well  within  the  typical  data  scatter. 

In  fact,  the  only  significant  observed  particle  velocity  which  is  not 
accounted  for  by  the  Rayleigh  wave  is  the  airslap-induced  motion  which  is 
identified  on  Figure  4  by  the  vertical  arrows  which  have  been  aligned  with 
the  predicted  airblast  arrival  times  at  these  ranges.  In  this  overpressure 
regime  (<  100  kPa),  the  ground  motion  associated  with  this  loading  can  be 
simulated  with  reasonable  accuracy  using  a  simple,  linear  one-dimensional 
model  which  incorporates  the  reverberation  of  energy  within  the  surface 
layer.  That  is,  if  p(t)  denotes  the  overhead  airslap  at  the  site,  then 
the  near-surface  particle  velocity  waveform,  U(t),  associated  with  rever¬ 
beration  of  the  airslap-induced  pulse  in  a  surficial  layer  of  thickness 
H  can  be  approximated  by  series: 
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Figure  3.  Fundamental  mode  Rayleigh  wave  site  response  function  for 
Frenchman  Flat  model. 
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Figure  4.  Comparison  of  observed  (top)  and  synthetic  (bottom)  particle  velocity  waveforms  for 
TUMBLER  I.  Vertical  arrows  denote  predicted  airblast  arrival  times. 
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The  relative  contributions  of  the  predicted  airslap-induced  and  Rayleigh 
wave  components  at  the  482  m  station  are  shown  in  Figure  5  where  they  are 
compared  with  the  observed  data.  It  can  be  concluded  on  the  basis  of  this 
comparison  that  the  linear  superposition  of  these  two  components  can  ac¬ 
count  for  virtually  all  the  significant  observed  motion  from  TUMBLER  I  at 
these  ranges. 

There  remains,  however,  one  minor,  puzzling  aspect  of  the  com¬ 
parisons  shown  in  Figure  4,  namely  that  the  synthetics  consistently  under¬ 
estimate  the  amplitude  of  the  observed  principal  velocity  pulses  by  a 
small  amount.  Although  this  discrepancy  is  well  within  the  normal  data 
scatter,  it  is  surprising  in  that  the  synthetics  have  been  computed  assum¬ 
ing  a  perfectly  elastic  response  and  this  should  provide  a  conservative 
upper  bound  to  the  observations,  at  least  on  the  average.  Thus,  the  ob¬ 
served  discrepancy  suggests  the  possibility  that  other,  unaccounted  for 
arrivals  may  be  making  noticeable  contributions  to  the  observed  motion. 

In  order  to  test  this  hypothesis,  the  contribution  of  the  first  higher 
mode  Rayleigh  wave  has  been  computed  and  compared  to  that  of  the  funda¬ 
mental.  Figure  6  shows  such  a  comparison  of  the  computed  fundamental  and 
first  higher  mode  contributions  for  the  407  m  station  together  with  the 
result  of  superposing  of  the  two  components.  It  can  be  seen  that  the 
higher  mode  waveform  is  more  impulsive  than  that  associated  with  the  funda 
mental  and  contributes  significantly  to  the  superposition  only  in  a  narrow 
arrival  time  window  centered  at  about  1  second.  Moreover,  comparing  the 
synthetic  superposed  waveform  of  Figure  6  with  the  waveform  observed  at 
that  range  from  Figure  4,  it  can  be  seen  that  the  first  higher  mode  contri 
bution  introduces  a  high  frequency  character  which  is  not  evident  in  the 
observed  signal.  However,  it  may  be  that  the  incorporation  of  some 
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Figure  5.  Comparison  of  observed  (top)  with  synthetic  (bottom)  airslap- 
induced  and  Rayleigh  wave  vertical  particle  velocity  waveforms 
for  TUMBLER  I,  R  =  482  m. 
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Comparison  of  fundamental  and  first  higher  mode  Rayleigh  wave 
contributions  to  vertical  particle  velocity  waveforms  for 
TUMBLER  I,  R  =  407  m. 


-'ealistic  anelastic  attenuation,  in  conjunction  with  the  addition  of  the 
higher  mode  contribution,  might  provide  a  somewhat  better  fit  to  the  obser 
vations.  In  any  case,  the  evidence  seems  to  indicate  that  the  higher  mode 
contributions  are  relatively  minor  in  this  case. 

2.2.2  UPSHOT-KNOTHOLE  10. 

1/3 

This  test,  with  its  scaled  HOB  of  about  65  m/kt  '  ,  is  inter¬ 
mediate  between  the  TUMBLER  I  and  the  near-surface  SMALL  BOY  explosions, 
and  is  unusual  in  that  it  was  fired  from  an  artillery  piece  (cf.  Table  1). 
Ground  motion  data  from  this  test  have  been  processed  and  published  by 
Stubbs  (1977)  and  the  vertical  particle  velocity  waveforms  recorded  at 
five  stations  in  the  horizontal  distance  range  extending  from  about  100  to 
1000  m  are  reproduced  in  Figure  7.  It  can  be  seen  that  the  data  are  of 
variable  quality,  but  that  all  exhibit  some  long  period  drift  which  makes 
it  difficult  to  resolve  and  analyze  the  lower  frequency,  oscillatory  com¬ 
ponents  of  the  recorded  ground  motion.  This  is  a  common  problem  with  much 
of  the  available  ground  motion  data  from  atmospheric  nuclear  explosions. 
Although  the  signals  of  Figure  7  have  been  digitized,  no  filtering  other 
than  a  simple  baseline  adjustment  has  been  applied,  so  as  to  preserve  the 
data  in  its  original  state.  While  this  is  a  commendable  approach  from  a 
documentation  standpoint,  it  seems  clear  that  the  utility  of  such  data  for 
purposes  of  ground  motion  analysis  could  be  greatly  increased  through  the 
application  of  suitable  highpass  filters.  Unfortunately,  the  digital  data 
corresponding  to  the  waveforms  in  Figure  7  are  not  readily  accessible  at 
the  present  time.  However,  for  purposes  of  illustration,  these  published 
waveforms  have  been  hand-digitized  and  filtered  using  a  simple,  second- 
order  highpass  filter  with  a  corner  frequency  of  0.5  Hz.  This  filter 
effectively  eliminates  spectral  components  with  periods  longer  than  about 
2  seconds,  thereby  improving  the  resolution  of  the  transient  ground  mo¬ 
tions.  This  is  illustrated  in  Figures  8-12  which  show  comparisons  of  the 
unfiltered  and  filtered  waveforms  for  the  five  stations  of  Figure  7.  As 
expected,  these  filtered  versions  provide  a  clearer  representation  of  the 
oscillatory,  low  frequency  motions  of  interest  in  the  present  investiga¬ 
tion.  Of  course,  these  results  need  to  be  treated  as  preliminary  in  that 
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the  resolution  which  can  be  obtained  by  hand-digitizing  waveforms  plotted 
to  the  scale  of  Figure  7  is  obviously  quite  limited.  This  is  particularly 
true  for  the  nearer  stations  where  the  motion  tends  to  be  dominated  by  the 
high-frequency,  airslap-induced  components.  Clearly,  much  more  reliable 
results  are  attainable  if  the  original  digitizations  of  Stubbs  (1977)  can 
be  recovered  and  filtered  in  this  manner.  In  the  meantime,  the  present 
simulation  analysis  of  UPSHOT-KNOTHOLE  10  will  focus  on  the  filtered  veloc 
ity  data  for  the  two  more  distant  stations  (i.e.,  R  =  736  m,  1041  m). 

As  with  TUMBLER  I,  the  low  frequency  ground  motions  from  UPSHOT 
KNOTHOLE  10  have  been  theoretically  simulated  using  the  Frenchman  Flat 
geologic  model  of  Figure  1  together  with  an  airblast  loading  predicted  by 
the  Speicher  and  Brode  (1981)  analytic  approximation.  Again,  as  a  first 
approximation,  rg  in  Equation  (1)  is  assumed  to  be  zero.  Figure  13  shows 
the  comparison  between  the  observed  and  fundamental  mode,  synthetic  verti¬ 
cal  particle  velocity  waveforms  for  the  most  distant  station  at  a  range  of 
1041  m.  It  can  be  seen  that  the  predicted  and  observed  waveforms  agree 
quite  well  with  respect  to  arrival  time,  dominant  frequency  and  amplitude 
level.  However,  the  onset  of  the  observed  Rayleigh  wave  arrival  is  more 
impulsive  than  that  of  the  synthetic,  suggesting  that  other  arrivals  may 
be  making  a  noticeable  contribution  in  this  case.  This  hypothesis  is  con¬ 
firmed  in  Figure  14  which  shows  a  comparison  of  the  predicted  contribu¬ 
tions  of  the  fundamental  and  first  higher  modes  at  the  range  of  1041  m. 

It  can  be  seen  that  although  the  1  Hz  oscillatory  component  of  the  first 
higher  mode  is  quite  small  relative  to  that  associated  with  the  fundamen¬ 
tal  mode,  the  initial  impulsive  motions  for  the  two  modes  are  of  compar¬ 
able  magnitude.  The  synthetic  obtained  by  superposing  the  contributions 
of  these  two  modes  is  shown  in  Figure  15,  where  it  is  compared  with  the 
observed  data.  Comparing  with  Figure  13,  it  can  be  seen  that  the  fit  to 
the  onset  of  the  surface  wave  arrival  has  been  significantly  improved  by 
the  addition  of  the  higher  mode  contribution.  In  fact,  as  with  TUMBLER  I, 
essentially  all  of  the  significant  motion  at  this  range,  with  the  excep¬ 
tion  of  the  airslap-induced  component  which  is  identified  by  the  vertical 
arrow  in  Figure  15,  can  be  explained  by  the  airblast-induced  Rayleigh  wave 
On  the  other  hand,  the  predicted  oscillatory  component  of  the  motion 
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Comparison  of  observed  and  synthetic  (fundamental  mode)  vertical  particle  velocity  waveforms 
for  UPSHOT-KNOTHOLE  10,  R  =  1041  m.  Vertical  arrow  denotes  predicted  airblast  arrival  time. 


appears  to  have  a  significantly  longer  duration  than  the  observed.  It  is 
not  clear  at  this  time  whether  this  discrepancy  is  due  to  problems  with 
the  observed  data  at  late  times  or  whether  it  is  a  real  effect,  associated 
perhaps  with  scattering  of  Rayleigh  waves  by  lateral  heterogeneities  along 
the  propagation  path.  A  careful  re-evaluation  of  the  original  digital  data 
may  help  to  resolve  this  issue. 

The  comparison  between  the  observed  and  synthetic  (fundamental 
plus  first  higher  mode)  vertical  particle  velocity  waveforms  at  the  736  m 
range  is  shown  in  Figure  16.  Again,  the  agreement  between  the  predicted 
and  observed  waveforms  in  terms  of  arrival  time,  frequency  content  and 
amplitude  level  is  generally  satisfactory.  However,  in  this  case,  the  ob¬ 
served  motion  is  dominated  by  a  secondary  high  frequency  arrival  which 
lags  the  predicted  airblast  arrival  time  (vertical  arrow)  by  more  than  1 
second  and  is  not  associated  with  the  Rayleigh  wave.  Careful  examination 
of  Figure  7  and,  more  particularly,  the  corresponding  acceleration  record¬ 
ings,  indicates  that  this  secondary  high  frequency  arrival  can  be  identi¬ 
fied  on  each  of  the  five  recordings.  The  source  of  this  arrival  has  not 
been  isolated  at  the  present  time,  but  one  possibility  is  that  it  is  asso¬ 
ciated  with  the  firing  of  the  280  mm  artillery  piece  that  was  used  to  fire 
the  explosive  to  the  desired  height  in  this  test.  This  hypothesis  can 
presumably  be  tested  if  the  location  of  the  gun  with  respect  to  the  record¬ 
ing  stations  can  be  recovered  from  the  test  documentation.  In  any  case, 
given  the  questionable  quality  of  the  available  data,  the  agreement  be¬ 
tween  the  observed  and  synthetic  particle  velocity  waveforms  at  these  dis¬ 
tances  from  the  UPSHOT-KNOTHOLE  10  test  is  quite  reasonable  and  confirms 
the  fact  that  much  of  the  observed  motion  can  be  accounted  for  by  the  air- 
blast-induced  Rayleigh  wave. 

2.2.3  SMALL  BOY. 

SMALL  BOY  was  a  low  yield  (i.e.,  <  20  kt)  explosion  detonated 
at  a  small  HOB  from  a  tower  over  Frenchman  Flat.  In  this  case,  it  seems 
clear  that  there  was  strong  nonlinear  interaction  between  the  airblast  and 
the  surface  area  around  ground  zero  and,  consequently,  a  nonzero  exclusion 
radius  (rn)  is  required  in  the  evaluation  of  Equation  (1).  Although  this 
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Figure  16.  Comparison  of  observed  and  synthetic  (fundamental  plus  first  highe 
mode)  vertical  particle  velocity  waveforms  for  UPSHOT-KNOTHOLE  10, 
R  =  736  m.  Vertical  arrow  denotes  the  predicted  airblast  arrival 


radius  has  been  found  to  coincide  approximately  with  the  final  crater 
radius  for  HE  surface  bursts  (Murphy,  1978),  there  is  no  available  evi¬ 
dence  to  suggest  how  it  should  vary  for  nuclear  explosions  with  nonzero 
HOB.  Consequently,  for  simulation  purposes,  a  preliminary  exclusion  radius 
of  30  m  around  ground  zero  has  been  somewhat  arbitrarily  adopted  to  define 
the  area  inside  of  which  no  airblast  loading  is  applied.  The  significance 
of  this  exclusion  radius  with  respect  to  low  frequency  ground  motion  sim¬ 
ulation  will  be  addressed  experimentally  later  in  this  section. 

Figure  17  shows  a  comparison  of  the  observed  SMALL  BOY  vertical 
displacement  at  a  range  of  2139  m  with  the  fundamental  mode  synthetic  com¬ 
puted  for  that  range  using  the  Frenchman  Flat  geologic  model  of  Figure  1. 

It  can  be  seen  that  the  agreement  is  not  very  satisfactory.  In  particu¬ 
lar,  although  the  amplitude  level  is  approximately  correct,  the  group 
arrival  time,  dominant  frequency  and  duration  of  motion  are  clearly  in¬ 
correct,  indicating  that  the  subsurface  model  of  Figure  1  is  not  an  ade¬ 
quate  representation  of  the  path  to  this  station.  In  fact,  it  is  difficult 
to  see  how  any  stronglv  layered  subsurface  model  could  give  rise  to  the 
observed  simple  displacement  pulse  at  this  relatively  large  range.  This 
observation  motivated  an  investigation  into  the  possibility  that  a  simple 
halfspace  model  might  provide  a  reasonable  approximation  to  the  observed 
data.  In  particular,  a  halfspace  model  was  selected  which  has  a  shear 
wave  velocity  consistent  with  the  observed  Rayleigh  wave  group  arrival 
time  at  this  station  (i.e.,  8  =  327  m/sec).  Figure  18  shows  a  comparison 
of  the  vertical  displacement  computed  using  this  halfspace  model  with  the 
observed  SMALL  BOY  data  at  a  range  of  2139  m.  It  can  be  seen  that  the 
agreement  between  theoretical  and  observed  obtained  using  this  model  is 
quite  good  with  respect  to  the  dominant  frequency,  amplitude  level  and 
duration  of  motion.  The  corresponding  particle  velocity  comparison  is 
shown  in  Figure  19.  Again,  the  agreement  is  very  good,  even  including 
the  relative  arrival  time  and  approximate  amplitudes  of  two  prominent, 
high  frequency,  secondary  arrivals  at  around  7.0  seconds.  It  is  somewhat 
puzzling  that  a  different  subsurface  model  is  required  to  explain  the  ob¬ 
served  ground  motion  from  SMALL  BOY,  given  its  proximity  to  the  other  two 
tests  (cf.  Table  1).  However,  it  was  detonated  near  the  center  of  the  dry 
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lake  bed  and  it  may  well  be  that  the  propagation  path  to  this  station 
crosses  a  different  portion  of  Frenchman  Flat  than  those  characteristic  of 
the  other  two  tests.  Unfortunately,  the  available  documents  describing 
the  ground  motion  experiments  on  these  tests  do  not  specify  the  map  coordi¬ 
nates  of  the  stations,  so  it  is  not  possible  to  compare  the  propagation 
paths  in  any  detail  at  this  time.  In  any  case,  the  observed  low  frequency 
ground  motions  clearly  indicate  that  such  differences  exist. 

Returning  now  to  the  particle  velocity  comparison  of  Figure  19, 
it  can  be  seen  that  the  theoretical  Rayleigh  waveform  is  considerably  more 
complex  than  that  which  would  be  expected  from  a  point  surface  load  on  a 
halfspace.  It  follows  that  the  propagating  airblast  loading  must  result 
in  a  complex,  multi -component  seismic  source  function.  For  an  observation 
point  at  large  distance  from  ground  zero,  the  Rayleigh  wave  source  function 
is  given  by  the  integral  expression  (Murphy  and  Bennett,  1980): 

/rmax 

p(r',w)  Jg(kr')  r'  dr'  (4) 

r0 

where,  for  the  applications  shown  in  Figures  18  and  19,  rQ  =  30  m  and 
r  v  =  1000  m.  Of  course,  by  linearity,  the  integral  in  Equation  (4)  can 

iTla  X 

be  decomposed  into  a  sum  of  integrals  of  the  form: 


(5) 


and  thus  it  is  possible  to  evaluate  the  contribution  of  any  arbitrary  sur¬ 
face  ring  loading  to  the  computed  waveform.  Figure  20  shows  a  comparison 
of  the  synthetic  particle  velocity  waveform  of  Figure  19  with  the  individ¬ 
ual  contributions  computed  from  the  loading  of  an  inner  ring  extending 
from  30  to  100  m  and  from  an  outer  ring  extending  from  100  to  1000  m. 

Note  that  the  two  high  frequency  secondary  arrivals  at  around  7.0  seconds 
arise  from  the  loading  of  the  inner  ring  (30-100  m),  while  the  initial 
lower  frequency  arrival  at  around  6.0  seconds  can  be  essentially  completely 
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accounted  for  by  the  contributions  from  the  outer  ring  (100-1000  m).  This 
is  physically  reasonable  in  that  in  the  inner  ring  the  overpressure  load¬ 
ing  is  impulsive  (i.e.,  has  a  short  duration)  and  highly  supersonic,  while 
in  the  outer  ring  the  positive  phase  duration  becomes  significant  (with 
respect  to  the  seismic  frequency  band  of  interest)  and  the  airblast  propa¬ 
gation  velocity  drops  down  into  the  seismic  velocity  range.  Another  ob¬ 
servation  that  can  be  made  with  regard  to  Figure  20  is  that  there  are  addi¬ 
tional  high  frequency  pulses  arriving  at  about  7.5  seconds  on  the  waveforms 
computed  from  the  individual  ring  loads  which  appear  to  have  no  counter¬ 
part  on  the  composite  waveform.  These  pulses  have  no  physical  significance 
and  are  merely  an  artifact  of  the  truncation  of  the  loading  at  100  m.  This 
sudden  truncation  introduces  a  negative  pulse  from  the  inner  ring  and  an 
equal  and  opposite  positive  pulse  from  the  outer  ring  which  cancel  out  in 
the  composite,  as  required  by  the  physics  of  the  problem. 

The  decomposition  of  the  loading  shown  in  Figure  20  also  pro¬ 
vides  insight  into  the  significance  of  the  exclusion  radius,  rg.  That  is, 
a  reduction  in  rQ  can  be  expected  to  increase  the  high  frequency  components 
of  the  induced  ground  motion  without  significantly  affecting  the  low  fre¬ 
quency  components.  This  inference  is  confirmed  in  Figure  21  which  shows 
a  comparison  of  the  observed  particle  velocity  waveforms  from  Figure  19 
with  the  synthetics  computed  assuming  rg  values  of  30  and  20  m  respectively. 
As  expected,  the  only  significant  change  in  the  synthetics  is  in  the  rela¬ 
tive  amplitude  of  the  high  frequency  secondary  arrivals.  In  fact,  this 
comparison  confirms  that  the  initial  value  of  30  m  is  about  right  and  that 
smaller  values  of  rg  lead  to  poorer  agreement  with  the  observations.  If 
anything,  a  slightly  larger  value  of  rg  is  suggested  by  these  data.  In 
fact,  nearly  exact  agreement  can  be  obtained  by  tapering  the  onset  of 
loading  at  30  m  rather  than  turning  it  on  discontinuously.  For  example. 
Figure  22  shows  a  comparison  of  the  observed  particle  velocity  data  with 
the  synthetic  computed  by  linearly  tapering  the  applied  overpressure  in 
the  30  to  60  m  range.  These  examples  confirm  that  the  exclusion  radius 
concept  is  a  reasonable  approximation  which  is  remarkably  consistent  with 
these  observations. 
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SECTION  3 


PREDICTION  AND  SCALING  CONSIDERATIONS 

3.1  BACKGROUND. 

The  results  presented  in  Section  II  of  this  report  have  con¬ 
firmed  the  fact  that  the  low  frequency  ground  motions  observed  from  se¬ 
lected  atmospheric  nuclear  tests  at  NTS  can  be  largely  accounted  for  by 
the  airblast-induced  elastic  surface  waves.  This  validation  under  full 
scale  nuclear  conditions  represents  a  significant  advance  over  the  pre¬ 
vious  analyses  of  HE  data.  However,  as  in  the  HE  studies,  the  simulation 
analyses  of  these  nuclear  events  were  "post-dictions"  of  data  which  had 
been  recorded  many  years  before  and,  thus,  not  directly  applicable  to  the 
assessment  of  prediction  capability.  Therefore,  in  the  next  section,  an 
example  of  some  recent  experience  with  the  actual  prediction  of  the  Pre- 
DIRECT  COURSE  experiment  low  frequency  ground  motions  will  be  presented 
to  illustrate  some  of  the  aspects  of  the  prediction  process.  This  will 
be  followed  by  an  outline  of  some  of  the  basic  scaling  relations  which  are 
implied  by  the  model  and  which  can  be  used  in  the  preliminary  evaluation 
of  the  effects  of  site  geology,  yield,  distance  and  instrument  depth  on 
the  predicted  low  frequency  components  of  the  ground  motion. 

3.2  Pre-DIRECT  COURSE  PREDICTIONS. 

This  experiment,  which  was  conducted  at  the  White  Sands  testing 
area  in  New  Mexico  in  September,  1982  consisted  of  24  tons  of  HE  (ANFO) 
detonated  at  an  HOB  of  about  17  m.  Preliminary  determinations  of  subsur¬ 
face  physical  properties  at  the  site  were  made  on  the  basis  of  a  P  wave 
refraction  survey  (Personal  Communication,  Bob  Reinke,  AFWL,  1982)  and  are 
shown  schematically  in  Figure  23  together  with  the  corresponding  predicted 
fundamental  mode  Rayleigh  wave  dispersion  curves.  It  should  be  noted  that 
the  shear  wave  velocities,  which  control  the  Rayleigh  wave  characteristics, 
have  not  been  measured  at  this  site  and,  consequently,  were  estimated  on 
the  basis  of  results  obtained  from  other  sites  in  this  general  area.  It 
can  be  seen  that  the  refraction  survey  indicates  the  presence  of  a  thick 
(75  m)  layer  of  dry  alluvium  overlying  saturated  alluvium  which  apparently 
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extends  to  a  depth  of  more  than  150  m  (i.e.,  the  approximate  limit  of 
resolution  of  the  500  m  long  refraction  line).  Note  that  because  of  the 
small  inferred  variation  of  shear  wave  velocity  (6)  with  depth,  very  strong 
dispersion  of  the  surface  wave  energy  is  not  expected  at  this  site.  How¬ 
ever,  the  model  does  predict  a  distinct  group  velocity  minimum  at  a  fre¬ 
quency  of  about  2.5  Hz  (i.e.,  8^/2H).  As  expected,  this  corresponds  to  a 
pronounced  corner  frequency  on  the  corresponding  Rayleigh  wave  site  response 
function,  Ar(ug),  shown  in  Figure  24.  Thus,  for  this  site  model,  the  pre¬ 
dominant  energy  in  the  surface  wave  components  of  the  ground  motion  should 
be  confined  to  frequencies  above  2.5  Hz. 

The  low  frequency  ground  motions  to  be  expected  from  this  test 
have  been  estimated  using  an  airblast  loading  function  computed  from  the 
Speicher  and  Brode  (1981)  analytic  approximation  for  nuclear  explosions, 
with  an  effective  yield  of  48  tons.  Available  evidence  suggests  that  a 
nuclear  explosion  of  this  yield  will  produce  an  airblast  loading  essen¬ 
tially  equivalent  to  that  expected  from  the  detonation  of  24  tons  of  ANFO 
at  that  same  height  (Personal  Communication,  George  Ullrich,  DNA,  1982). 
However,  in  the  case  of  Pre-DIRECT  COURSE  this  clearly  provides  only  a 
rough  approximation  to  the  loading  in  that  the  mass  of  the  container  was 
about  12  percent  of  the  total  explosion  mass  and  this  significantly  per¬ 
turbed  the  airblast  characteristics  in  the  high  overpressure  regime. 

Figure  25  shows  a  comparison  of  the  observed  Pre-DIRECT  COURSE 
vertical  displacements  at  ranges  of  94,  156  and  313  m  with  the  displace¬ 
ments  predicted  theoretically  using  the  geologic  model  of  Figure  23  to¬ 
gether  with  the  airblast  loading  approximation  described  above.  It  can  be 
seen  that  while  the  predicted  and  observed  amplitude  levels  agree  quite 
well,  the  predicted  dominant  frequency  is  too  low  and  the  predicted  arri¬ 
val  time  is  consistently  late.  This  suggests  that  the  geologic  model  of 
Figure  23  is  not  an  accurate  representation  of  the  actual  subsurface  con¬ 
ditions  at  this  site.  Consequently,  the  approximate  scaling  law  f  = 

3j/2H  (Murphy  and  Bennett,  1980)  has  been  used  to  infer  a  model  which  is 
more  consistent  with  the  observed  ground  motion  characteristics.  The 
final  model  is  compared  with  the  original  model  in  Figure  26  where  it  can 
be  seen  that  the  principal  change  has  been  to  decrease  the  average 
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surficial  layer  thickness  from  75  m  to  50  m.  The  theoretical  displace¬ 
ments  predicted  using  this  final  geologic  model  are  compared  with  the  ob¬ 
served  data  in  Figure  27.  It  can  be  seen  that  the  agreement  is  now  well 
within  the  variability  in  the  experimental  data.  Moreover,  subsequent  de¬ 
tailed  geophysical  surveys  at  the  site  revealed  that  the  actual  surficial 
layer  thickness  is  indeed  about  50  m  as  opposed  to  the  original  estimate 
of  75  m  (Personal  Communication,  George  Ullrich,  DNA,  1983).  Thus,  the 
evidence  provided  by  this  experiment  suggests  'hat  given  an  accurate  model 
of  the  subsurface  site  geology,  it  should  be  possible  to  accurately  pre¬ 
dict  the  characteristics  of  the  low  frequency  (round  motions  expected  from 
atmospheric  explosions.  However,  the  predicteu  gr'»u'-u  motion  characteris¬ 
tics  are  clearly  quite  sensitive  to  the  near-surface  shear  wave  velocity 
structure,  and  the  Pre-DIRECT  COURSE  experience  indicates  that  this  cannot 
be  adequately  constrained  on  the  basis  of  P  wave  refraction  surveys. 

Another  conclusion  that  can  be  drawn  from  this  analysis  is  that  departures 
of  the  close-in  airblast  loading  from  the  ideal  represented  by  the 
Speicher/Brode  analytic  approximation,  such  as  those  which  occurred  on 
Pre-DIRECT  COURSE,  do  not  have  a  pronounced  effect  on  the  observed  low 
frequency  ground  motions.  That  is,  as  has  been  noted  previously  (Murphy 
et  al_.,  1982),  these  ground  motion  components  are  relatively  insensitive 
to  the  high  frequency  details  of  the  airblast. 

3.3  SCALING  CONSIDERATIONS. 

The  results  presented  above  and  in  the  previous  section  indi¬ 
cate  that  it  is  possible  to  accurately  estimate  the  complete  waveforms 
associated  with  the  low  frequency  components  of  the  ground  motions  induced 
by  atmospheric  explosions,  given  the  yield  and  HOB  of  the  explosion  and 
an  accurate  description  of  the  subsurface  site  geology.  However,  in  many 
applications  in  which  the  ground  motion  vulnerabilities  of  various  strate¬ 
gic  structures  are  being  parametrically  evaluated  in  a  preliminary  fashion, 
such  detailed  predictions  are  not  required  and  what  is  needed  are  simpli¬ 
fied  scaling  laws,  such  as  those  provided  in  EM-1.  At  the  present  time, 
work  is  still  in  progress  on  the  development  of  a  simple,  theoretically- 
based  scaling  model  which  can  be  used  to  approximate  all  the  low  frequency 


ground  motion  parameters  of  interest.  However,  the  basic  scaling  laws  de¬ 
scribing  the  variation  of  peak  amplitude  and  dominant  frequency  with  yield, 
range,  site  geology  and  instrument  depth  have  already  been  inferred  from 
the  model  and  provide  some  preliminary  basis  for  parametric  analysis. 

The  principal  assumptions  underlying  the  theoretical  scaling 
laws  may  be  enumerated  as  follows: 

(1)  The  subsurface  geology  at  the  site  of  interest  can 
be  approximated  as  a  sequence  of  homogeneous  layers 
of  constant  thickness. 

(2)  The  source  of  the  surface  waves  can  be  described  by 
the  linear  response  of  the  medium  to  the  airblast 
loading  on  the  surface  outside  the  crater  region 
(if  any) . 

(3)  The  yields  of  interest  are  sufficiently  large  that 
the  surface  wave  motions  will  be  dominated  by  fre¬ 
quency  components  near  the  corner  frequency  of  the 
Rayleigh  wave  site  response  function,  which  coincides 
with  the  frequency  of  minimum  group  velocity  and  is 
given  approximately  Dy  wnere  h  denotes  tne 
total  soil  thickness  at  the  site  and  B  is  the  aver¬ 
age  shear  wave  velocity  of  the  soil  column  (Murphy 
and  Bennett,  1980). 

(4)  The  observation  point  of  interest  is  close  enough 
to  the  source  that  the  effects  of  anelastic  attenu¬ 
ation  on  the  propagating  surface  wave  can  be 
ignored  to  first  order. 

Given  these  assumptions,  the  surface  wave  peak  amplitude.  A,  at  a  parti¬ 
cular  site  is  predicted  to  scale  in  the  following  manner. 

•  At  a  fixed  observation  distance  (Murphy  and  Bennett,  1980): 


At  a  fixed  yield,  the  amplitude  decay  with  distance, 
r,  is  given  approximately  by  the  Airy  phase  attenua¬ 
tion  law: 

A  -  r‘5/6 

At  a  fixed  scaled  distance,  it  follows  from  the  above 


two  laws  that: 


A  ~  w 


Note  that  under  the  assumption  that  the  surface  wave  is  dominated  by  a 
pure  sinusoidal  component  with  constant  frequency  B/2H,  these  scaling  laws 
apply  to  the  velocity  and  acceleration,  as  well  as  the  displacement,  peak 
ampl itudes. 

The  expected  depth  dependence  of  the  surface  wave  motions  can 
be  analyzed  in  a  preliminary  fashion  by  employing  the  same  approximations 
used  in  defining  the  other  scaling  laws.  Let  Uj(Z)  and  W^Z)  denote  the 
horizontal  and  vertical  Rayleigh  wave  displacement  components  as  a  func¬ 
tion  of  depth  in  the  "surface  layer".  Then  for  cases  in  which  the  shear 
wave  velocity  ratio  <  4  we  can  write  for  frequency  components  near 

3-.  /2H  (Murphy  and  Hewlett,  1975): 
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and  c  is  the  phase  velocity  at  angular  frequency  u.  For  the  nominal  case 
in  which  a  =  /J  3  and  3^  ~  23^,  it  can  be  shown  that  at  the  Rayleigh  wave 
characteristic  frequency  a>  ~  ttBj/H,  c/31  »  1.3,  £  =  0.6  and  Equations  (6) 
and  (/)  can  be  rewritten  explicitly  in  tne  torm: 


Wj(Z) 

Wj(0) 


=  0.48  sinh  1.43  J  -  0.21  cosh  1 .43  £ 
-  0.15  sin  1.72  J  +  1.21  cos  1.72  jj 


U! ^ z)  z  z 

— -  =-0.71  cosh  1.43  n  +  0.31  sinh  1 .43  £ 

iWj(0)  H  H  (14 

+  0.12  cos  1.72  jj  +  0.98  sin  1.72  \ 

The  depth  dependences  predicted  by  these  equations  are  shown  in  Figure  28 
as  a  function  of  the  dimensionless  variable  Z/H.  It  can  be  seen  that,  as 


in  the  case  of  the  halfspace,  the  horizontal  component  of  motion  decreases 
to  zero  fairly  rapidly  (i.e.,  at  Z  ^  0.35H)  and  then  changes  sign.  As  is 
indicated  on  the  figure,  these  depth  dependences  can  be  approximated  quite 
closely  by  the  following  simple  expressions. 
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where  ( U, (0) |  =  0. 6 | Wx (0) | . 

Preliminary  indications  are  that  the  theoretically-based  scaling 
laws  proposed  above  are  generally  consistent  with  the  large  body  of  low 
frequency  ground  motion  data  which  have  been  observed  from  both  HE  and 
nuclear  atmospheric  explosions  (Personal  Communication,  C.  J.  Higgins, 

ARA,  1984).  Work  is  currently  in  progress  to  perform  more  definitive  com¬ 
parisons  with  the  available  observations  and  to  define  the  range  of  appli¬ 
cability  of  these  approximate  relations. 


SECTION  4 


CONCLUSIONS 

4.1  SUMMARY. 

The  investigations  summarized  in  this  report  have  centered  on 
the  simulation  analysis  of  low  frequency  ground  motions  observed  from 
atmospheric  explosions  and  on  the  preliminary  evaluation  of  a  theoretically 
based  prediction  methodology.  In  particular,  the  present  effort  has  ex¬ 
tended  the  previous  research  on  low  frequency  ground  motions  from  HE  experi 
ments  to  a  full  scale  validation  study  conducted  using  data  recorded  from 
atmospheric  nuclear  explosions  at  NTS. 

The  mathematical  model  being  used  to  simulate  the  low  frequency 
ground  motions  produced  by  atmospheric  explosions  was  briefly  summarized 
in  Section  II  and  then  applied  to  the  theoretical  simulation  of  ground 
motion  data  recorded  from  the  NTS  nuclear  explosions  TUMBLER  I,  UPSHOT- 
KNOTHOLE  10  and  SMALL  BOY.  These  three  explosions  represent  a  range  of 
scaled  HOB  extending  from  about  240  m/kt^3  (TUMBLER  I)  to  essentially 
zero  (SMALL  BOY)  and  thus  provide  an  opportunity  to  assess  simulation 
capability  over  essentially  the  entire  range  of  interest  in  ground  motion 
vulnerability  studies.  The  results  of  these  simulations  have  confirmed 
the  fact  that  the  observed  low  frequency  ground  motions  from  such  explo¬ 
sions  can  be  accounted  for  by  the  airblast-induced  Rayleigh  wave.  More 
specifically,  these  analyses  have  provided  additional  insight  into  the 
variation  of  the  Rayleigh  wave  seismic  source  function  with  HOB.  Thus, 
the  simulation  analysis  of  the  ground  motion  data  recorded  from  the  near¬ 
surface  SMALL  BOY  explosion  confirmed  the  existence  of  an  "exclusion  radius 
for  such  events,  within  which  the  incident  airblast  is  not  efficiently 
coupled  into  the  Rayleigh  wave  radiation  field  due  to  energy  dissipation 
associated  with  strong,  nonlinear  interaction  effects.  On  the  other  hand, 
the  evidence  provided  by  the  TUMBLER  I  and  UPSHOT-KNOTHOLE  10  analyses  in¬ 
dicated  that  this  exclusion  radius  is  effectively  zero  for  explosions  with 
significant  HOB. 

Some  preliminary  considerations  concerning  the  prediction  of 
low  frequency  ground  motion  parameters  were  addressed  in  Section  III, 
where  the  experience  gained  with  the  Pre-DIRECT  COURSE  experiment  was 


49 


assessed  and  some  approximate  scaling  laws  describing  the  dependence  of 
low  frequency  ground  motion  characteristics  on  variables  such  as  yield, 
distance,  site  geology  and  instrument  depth  were  presented  and  discussed. 

It  was  found  that  the  initial  pre-shot  prediction  of  the  Pre-DIRECT  COURSE 
low  frequency  ground  motions  was  only  partially  successful,  due  to  the 
fact  that  the  subsurface  geologic  model  available  at  that  time  was  inaccu¬ 
rate.  Subsequent  analyses  of  the  measured  data,  using  a  modified  site 
model  inferred  from  more  extensive  post-shot  exploration,  provided  evidence 
which  indicated  that  it  should  be  possible  to  predict  the  low  frequency 
ground  motions  produced  by  atmospheric  explosions  with  good  accuracy  using 
the  proposed  theoretical  model,  if  an  accurate  subsurface  site  model  can 
be  specified. 

4.2  CONCLUSIONS. 

The  analyses  summarized  above  support  the  following  conclusions 
regarding  the  low  frequency  ground  motions  observed  from  atmospheric  nu¬ 
clear  explosions  at  NTS. 

(1)  The  low  frequency  components  of  the  ground  motions 
observed  from  these  explosions  can  be  accounted  for 
by  the  airblast-induced  Rayleigh  waves.  In  fact, 
in  the  regime  where  the  overpressure  is  less  than 
about  100  kPa,  this  component  of  the  motion  was 
found  to  account  for  essentially  all  the  significant 
observed  particle  velocity  excursions  with  the  ex¬ 
ception  of  that  associated  with  the  overhead  airslap- 
induced  motion. 

(2)  There  is  an  "exclusion  radius"  associated  with  near¬ 
surface  explosions,  inside  of  which  incident  energy 
is  not  efficiently  coupled  into  the  surface  wave 
mode  of  propagation.  This  radius  has  been  shown  to 
be  effectively  zero  for  explosions  with  significant 
HOB. 


(3)  Much  valuable  ground  motion  data  have  been  recorded 
from  atmospheric  nuclear  explosions.  However,  these 
data  need  to  be  better  documented  and  reprocessed 
using  modern  signal  analysis  techniques  so  that  they 
can  be  more  fully  utilized  in  ground  motion  analysis 
studies . 

With  regard  to  development  of  procedures  for  predicting  the 
characteristics  of  the  low  frequency  ground  motions  produced  by  atmospheri 
explosions,  results  of  studies  conducted  to  date  lead  to  the  following 
concl us  ions . 

(1)  The  predicted  surface  wave  characteristics  are  quite 
sensitive  to  the  near-surface  shear  wave  velocity 
structure  at  the  site.  Pre-DIRECT  COURSE  experience 
indicates  that  this  velocity  distribution  cannot  be 
adequately  constrained  on  the  basis  of  P  wave  refrac¬ 
tion  surveys. 

(2)  For  explosions  which  are  large  enough  that  the  sur¬ 
face  wave  motions  are  dominated  by  frequency  compo¬ 
nents  near  the  corner  frequency  of  the  Rayleigh  wave 
site  response  function  (i.e.,  f  =  S/2H),  the  theo¬ 
retical  model  predicts  that  the  amplitude  of  the 

low  frequency  ground  motion  should  vary  with  yield, 

1  /3 

W,  as  W  and  with  range,  R,  approximately  as 
R-5'6. 

(3)  If  W(0),  U(0)  denote  the  vertical  and  horizontal 
component  Rayleigh  wave  displacements  on  the  surface, 
characterized  by  dominant  frequency  3/2H,  then  the 
corresponding  displacements  at  depth  Z  within  the 
surface  layer  of  total  thickness  H  can  be  approxi¬ 
mated  by  the  relations: 


SECTION  5 


Hawthorne 

Murphy,  J 

Murphy,  J 

Murphy,  J 

Murphy,  J 

Murphy,  J 
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